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A MATRIX MODEL FOR SIMPLE HURWITZ NUMBERS, AND 
TOPOLOGICAL RECURSION 

GAETAN BOROT, BERTRAND EYNARD, MOTOHICO MULASE, AND BRAD SAFNUK 

Abstract. We introduce a new matrix model representation for the gener- 
ating function of simple Hurwitz numbers. We calculate the spectral curve 
of the model and the associated symplectic invariants developed in [4]. As 
an application, we prove the conjecture proposed by Bouchard and Mariiio 
[2], relating Hurwitz numbers to the spectral invariants of the Lambert curve 



1. Summary 

In [2], Bouchard and Marino propose a new conjectural recursion formula to 
compute simple Hurwitz numbers, i.e. the weighted count of coverings of CP^ with 
specified branching data. Their recursion is based on a new conjectured formal- 
ism for the type B topological string on mirrors of toric Calabi-Yau threefolds, 
called "remodeling the B-model", or "bkmp conjecture" [T]. The Bouchard-Marifio 
conjecture for Hurwitz numbers appears as a consequence of this general bkmp 
conjecture applied to the infinite framing limit of the open string theory of C'^. In 
this limit, the amplitudes are known to give simple Hurwitz numbers. 

They propose that the generating function for Hurwitz numbers can be recovered 
from the symplectic invariants (also called topological recursion) developed in [3], 
applied to the so called "Lambert curve" y = L{e^) defined by: 

In this paper, we make the link between Hurwitz numbers and the Lambert curve 
explicit. We introduce a new matrix model formula for the generating function of 
simple Hurwitz numbers 

(1) Z cx / dMexp ( - — Tr(F(M) - MRU , 

JuNiC) \ 9s J 

where V{x) is the potential 

1 

^{^) = -y +gs{N - -)x + x\n{gs/t) + iirx - ln(r(-a;/.gs)) . 

The parameters gs and the matrix R involved in the definition of Z are such 
that the weight of a covering of Euler characteristic x is proportional to gf, and 
has a polynomial dependance in Vi = exp Ri which encodes the ramification data 
above branch points. 

A method to compute topological expansion of matrix integrals with an external 
field was introduced in [4]. It consists in finding the spectral curve S (roughly 
speaking the equilibrium density of eigenvalues of the matrix, more precisely the 
planar part of the expectation value of the resolvent), then computing recursively a 
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sequence of algebraic /c-forms yvjf\s), and some related algebraic quantities called 



symplectic invariants Tg{S) = wlf\s). Then, one of the main results of ,4, is that 



a=o 

In our case, this implies that the generating function for simple Hurwitz numbers 
of genus g is precisely !Fg (S) , where S is the spectral curve of our matrix model. 

It is rather easy to find the spectral curve of the matrix model Eqn[T] The result, 
after suitable symplectic transformations, reads 



(2) 5(p,<7,;t) 




z + In {z/t) + CO + f - EZi f2n{z) 



9sT, 



^ 1 I _J_ 

= 1 {z-Zi)yi ~^ ZiVi 



where Zi, yi, Cq and ci are determined by consistency relations, and Zi and are 
0(1) when — > 0, and Cq and Ci are 0{gs). In particular, when we set the coupling 
constant gs = 0, we recover the Lambert curve y = L{te^). 

In [2], Bouchard and Marino define another set of generating functions, denoted 
H^^'^{xi, . . . , Xk), encoding genus g simple Hurwitz numbers, and which are deriva- 
tives of In Z, evaluated at gs — 0. The statement of their conjecture is: 

dx{zi) ■ ■■dx{zk) 

where yvjf^ are the k-forms of U computed for the Lambert curve. 

We prove their conjecture by using the general properties of the invariants of 
[Ij, in particular the fact that derivatives of the Tg's with respect to almost any 

parameter, can be expressed in terms of the wife's. Then it suffices to set g^ = 0, 
and this gives the Bouchard-Marino conjecture. 

Organization of the paper. In Section |2] we recall the definitions and derive 
a matrix model formula for the generating function of simple Hurwitz numbers. 
We recall the construction of the symplectic invariants and topological recursion 
of [3] in Section [s] In Section |4] we derive the spectral curve of our matrix model 
(the proof is presented in Appendix [?]) and prove the Bouchard-Marino conjecture 
following a method very close to 0. In section [sj we briefly study the link with the 
Kontsevich integral. Section [6] addresses generalizations of our method and open 
questions. 



2. Construction of the matrix model 

Let Gov* (Ci, . . . , Cfe) denote the weighted number of n-fold coverings (possibly 
disconnected) of CP^ ramified over k fixed points of CP^ with monodromies in the 
conjugacy classes Ci , . . . , . The weight is one over the order of the automorphism 
group of the covering. Similarly, let Cov„(Ci, . . . , Ck) denote the weighted number 
of n-fold connected coverings. 

By a result of Burnside (see, eg [H]), we have 
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where the sum ranges over all partitions A = (Ai > A2 > • • • > A„ > 0) of 
|A| = = n boxes, dim A is the dimension of the irreducible representation 

indexed by A (with corresponding character xa), and 

/a(Q) = j^XA(C,). 
dim A 



in the 

n-fold covering 



in CP, 





Example of 
ramification profile 
above one point 
of monodromy class 



Example of a simple branch point 



Figure 1. Branched covering, with one branch point of monodromy class 
fi, and simple branch points (for which the monodromy is a transposition). 
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2.1. Simple Hurwitz numbers. We are interested in counting rt-fold coverings 
of genus g with 1 branch point of arbitrary profile /i, and only transpositions above 
other points (called simple branch points), c.f. Figure |lj We denote C(2) the 
conjugacy class of a transposition. For h simple branch points and one branch 
point of profile /i, the Euler characteristic of the n— fold covering reads from the 
Riemann-Hurwitz formula: 

x^\^Ji\ + ^{^i)-h. 

For connected coverings, we have x = 2 — 2g, where g is the genus, and we also 
define the simple Hurwitz numbers: 

6 

. ^ . 

Hg^p, = Cov„(C^, C(2), . . . , (^(2)), 

where & 2.g - 2 + |^| + 

2.2. Generating function for simple Hurwitz numbers. With the notations 
~ Y\iP ij.il E^J^d p — {pi,P2i ■ ■ ■), we shall study the generating function 

oo oo b-\^l\-e{^l) ^ 
(4) Z(p,5,;i) = ^r J2 T.-^i p^Cov:(C^,C(2),...,C(2)), 

n=0 \fJ.\=n b=0 

where Gov* (C^, C(2), . . . , (^(2)) is the number of (not necessarily connected) branched 
coverings of Euler characteristic X = l/^l + ^(m) ~ ^• 

In the language of string theory, —gs is the string couphng constaniQ Let us 
emphasize that Z{p,gs;t) is defined as a formal power series in t and g^, i.e. it is 
merely a notation to collect all the coefficients. Each coefficient (for b and n fixed) 
is a finite sum, which is a polynomial function of pi, . . . ,p„. Notice also that the 
parameter t is redundant because we can change pj — > p'pj and t t/p without 
changing the sum, i.e. 

Z{{pi,P2,P3, ■ ■ ■ ,},.9s;i) = Z{{ppi,p^p2,p^P3, . . .,},gs;t/p). 

In [5], t is chosen as t = 1, but we find more convenient to keep t 1 for the moment, 
in order to have only two formal parameters gs and t, instead of an infinite number 
of them gs and pi,p2, . . . , which would be the case if t were set to 1 . The generating 
function of connected coverings is F = In Z (in the sense of formal power series 
of t and gs): 

F(p, gs;t)=lnZ^Y.u ^ 

b.n \tJ\=n 

where = 2^ - 2 + |^| + ({p). 

Therefore, we have a so-called topological expansion (equality of formal series): 

00 

F{p,gs;t)=Y.9's''' ^^(P;*): 

3=0 

where Fg counts the number of connected coverings of genus g: 
Fgip; t) = E r ^£ (2,_2/: + £(^)). 



In the physics htcrature, it is customary to choose —gs instead of Qs as formal parameter. 



A MATRIX MODEL FOR SIMPLE HURWITZ NUMBERS, AND TOPOLOGICAL RECURSIONS 



Our goal in this article is to provide a recursive algorithm to compute the i^g's, 
and more precisely, prove that the Fg^s are the symplectic invariants introduced 
in |3] for a spectral curve S{p,gs;t) which we shall describe in Section [i] As a 
consequence, we shall prove the conjecture of Bouchard and Marino 



2.3. Partitions. To make our notations clear, we recall some representations of 
partitions or Young tableaux. The set of all partitions A of length < TV is in 
bijection with other interesting sets of objects: 

• The decreasing finite series of N integers : Ai > . . . > A^^;^) > Xi(^x)+i — 
... — Xn > 0. The length £{X) of the partition is the number of non- 
vanishing Xi's. |A| = '^i=i Ai = n is the number of boxes of the partition. 

• The strictly decreasing finite series of positive integers. They mark the 
positions (up to a translation) on the horizontal axis of the increasing jumps 
when the Young tableau {Xi,i) is tilted anticlockwise by j (c.f. Figure |2]). 
The correspondence is given by 

h, = K-i + N {ie{i,...,N}). 
We have /ii > /12 > ■ ■ • > > 0. 




Figure 2. Rotate the partition by 7r/4. Tlie h^'s mark the positions (up to 
a translation) on the horizontal axis of the increasing jumps in the Young 
tableau. 



• The conjugacy classes of 6„. The class C\ associated to A is the one with 
rrir = |{z > Xi — r}\ cycles of length r, and its cardinal is 

n, >i '"r! • r-™- ' 

• The equivalence classes of irreducible representations of ©|a|- 

• The equivalence classes of irreducible representations of GL((x){C), or of 
UWA)). 
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2.4. Schur polynomials. Recall the definition of Schur polynomials s\ [13j they 
coincide with the characters of the representation of U(fc) indexed by A such that 
k = £{X). As a matter of fact, if v = (vi, . . . , Vk) is an fc-tuple of complex variables: 



(5) 5a(v) = 



det («^-^+^ 

A(vr" 



where A(v) = Y[i<j<i<Ni''^i ~ ^j) Vandermonde determinant. This formula 

can be extended to a definition of sx with N > k variables, v = (wi, . . . , wat), 
provided that we take Xj = whenever j > k. 

The Frobenius formula gives the expansion of a Schur polynomial in terms of 
the power-sum functions pm — X^i^i ^I"; 

(6) sa(v) = ^ 5] |C^|xa(C^)p,„ 

which stresses the link between U(fc) characters and 6„ characters. It is still valid 
with N > k variables instead of ^(A) — k variables. 

From Eqn[5l one can obtain the classical result for the dimension of the repre- 
sentation indexed by A: 

sa(1, . . . , 1) = dim A = 



while the Frobenius formula leads to the expression 



2.5. Z as a sum on partitions. After Eqn|6] if we consider the variables Pm's to 
be power-sum functions of some iV-upple parameter v, we have 

(7) I ^^P'9s;t) = Eh,„ * 6! J2£{\)<N,\\\=n ■Sa(v) (/a(C(2)))^ 



V" 



Hence, 



Zip,g,;t)^ Y: ^ ^a(v) e^=^^(^<^)). 

i{X)<N ' 

Alternatively, given that s\ is homogeneous of degree |A|, we have 

i{\)<N ' 

Again, we emphasize that the above expression for Z ought to be considered as 
a formal power series in t and gs- For a given b and n, the coefficient of is 
a polynomial in the Pm's, which involves only pi, . . . ,p„. Therefore, it is always 
possible to find N (independent of b) and Wi, . . . , uat such that Vi S {!,..., N} we 
have Pm = gsJ2''^T- The values of the Pm's for m > n do not matter. 
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2.6. Matrix integral representation of Z. To write Z as a matrix integral we 
express s\ in terms of the Itzykson-Zuber integral |12j : 

_ det(e^'W) 

A(x)A(r)- 

Here dU is the Haar measure on U (N) , normalized such that the second line holds 
without any constant prefactor. 

Therefore, if we let R — diag(lni;i, . . . ,lnt;jv) and = diag(/ii, . . . , hj\j), we 
have from Eqn |5] : 

.A(v) = A(h,)|^/(h,,R). 



Then, the partition function looks like 

^ ' A i=l 



- ^ E nh.R)(A(h))'n ' r, f , — 

- MA(vy ^ /(h,R)(A(h)) _[| Til^) ' 

where |A| = X)i^i(^i) f\{C(2)) — J2i^^i^i)- The factorization of the weight 
with respect to the hi's and the presence of a squared Vandermonde is the key for 
the representation of Z as a hermitian matrix model. 

As it was done in [B], we represent the iV sums as integrals over a contour, 
namely the contour Cq enclosing the non-negative integers, as pictured in Figure [3] 
We make use of a function which has simple poles with residue 1 at all integers: 

/(O - - -ra + i)r(-Oe-"«. 

sni(7r4) 



Thus, we have: 



Z{p,gs;t) = 1^^^^ £jh,---dhN{A{h)yi{h,R) 

T{h, + 1) 



iV! A(v) 



Actually, the ratio /(^) /r(^ + 1) has only poles at non- negative integers, it has no 
pole at negative integers. So, we can replace Co by the contour C from Figure [3] 
which encloses all integers, and which we choose invariant under translation on the 
real axis. Therefore, we arrive at 



Z{p,gs;t) = ^^^^^d/ii---d/i^(A(h))'/(h,R) 

N 

n 



T{h, + 1) 
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CPi 



■< 
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■ * 





Figure 3. Contours enclosing the non- negative integers and the 
entire real line. 



The set Tipf^C) of normal matrices with eigenvalues on C is the set of matrices 
M which can be diagonalized by conjugation with a unitary matrix, and whose 
eigenvalues belong to C: 

M = U^XU, UU^ ^Id, X = diag(a;i,...,XAr), a:, G C. 

TLn{C) is endowed with the measure 

dM = A(X)2 dXdU, 

where dU is the (up to normalization) Haar measure on U (N) and dX is the product 
of Lebesgue curvilinear measures along C. 

From the above discussion, we can express our generating function as the matrix 
integral 

1 A(R) 



iV! A(v) Jh^{c) 



where 



v{0 = -gsA^iO + H9s/t) Ai(0 + - In (r(-e)). 

Since we are interested in the expansion as a power series in , we prefer to rescale 
^ = x/Qs, i.e. M M/gs, and rewrite 



(8) 



Zip,gs;t) 



A(R) 



A(v) y^„(c) 



dMe ? 



i Tr(v'(M)-MR) 



where now the potential reads 

= ~g'^^A2{x/gs) + gs ln{gs/t)Ai{x/gs) + inx - g^ In {r{-x/gs)) 

^-\+9s{N -]^)x + (Hgjt) + n:)x - g^ In {T{-x/gs)) + Ct, 
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with 



C, 



-Ig'siN' - + 2) + \g.{N - 1) ln(5.A) 



To write its derivative, let ip = T'/T, then 



V'{x) = -x + g,{N--)+ \n{gjt) +m + i){-x/g,). 



We have the well known formula that: 



where 7 is the Eulcr-Maschcroni constant. In other words, V'(0 "^'^^ simple polos at 
all negative integers, and an essential singularity (log singularity) at ^ = 00. This 
gives: 

1 1 °° / 1 1 \ 

(9) V'[x) = -x + g.iN--^) + Hgjt) + + 7 - ^ - g + ^ j • 

V'{x) has simple poles at x = Xj = jg^ for all j G N, and an essential singularity 

at a; = 00. 

Using Stirling's formula for the large ^ asymptotic expansion of tp 



2^ 



l> 



- 21 e 



we have order by order in the small gs expansion 

(10) V'{x) = -x + ln {x/t) +g^iN-^) + ^-J2^%- 

1=1 



where the -B^'s are the Bernoulli number. Wc recall their generating function 

e« - 1 ^^2^ n^- 



The first few are B2 = 1/6, B4 = -1/30, . . . while B21+1 = for / > 1. 



3. Generalities on matrix models 



Matrix integrals of the form 

Z = [ dMe-il^[^(^'-H 

are called " 1-matrix model in an external field" . They can be computed for any 
potential V, any contour C and any external matrix R. Here, our task is even 
simpler because we regard this integral as a formal integral, i.e. a formal power 
series in powers of t and gs, and therefore all our computations are to be performed 
order by order in powers of t and gs- 

In our case V depends on g,,, but for the moment, let us assume that V is 
an arbitrary potential, and in particular we assume that there is no relationship 
between the coefficients of V and gg. 
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Typically, in our case (see Eqn [9| , we choose V of the form: 



(11) V'{x) = ~x + C + Y,^^_, 



where for the moment we assume that there is no relationship between the coeffi- 
cients of Uj , Xj of V and gs ■ 

3.1. Topological expansion. For some choices of V, R and C, it may or may not 

happen that the convergent integral Z has a power series expansion in gs 

oo 

In Z = J2gl^-'F,. 

3=0 

This happens only if C is a "steepest descent path" for the potential V and R. 
In general, it is rather difficult to compute the steepest descent paths of a given 
arbitrary potential. 

Fortunately, when Z is defined as a formal integral we don't need to find the 
steepest descent paths, and very often, formal series do have a topological expansion 
almost by definition, order by order in the formal parameters, which is the case here 
as we argued in Paragraph |2.2[ In other words. 



lnZ(p,g,;t)=^5,2s-2^,(p;t) 

5=0 

holds order by order in powers of gs and t. 

3.2. Loop equations and spectral curve. We introduce the resolvent Wi and 
auxiliary quantities Pi and Pij'. 



Wi{x) = (Tr-^) = ^g,2^-iW^(«)(a;) 



- 3=0 



oo 

= (((x,-M)-%)=^5^^-i/^f 

9=0 



which we assume to have topological expansions in powers of g1^~^. In our case, 
their precise definition, order by order in i as a power series of gs, is given in 
Appendix [7] 

Loop equations, also called Schwinger-Dyson equations, is a general technique, 
which merely reflects the fact that an integral is invariant by change of variable. It 
is a standard matrix model exercise, (see ^ for the 1-matrix model in an external 
fiels), to prove that the invariance of Z under the infinitesimal change of variable 

x — M y ~ H, 
implies that the following loop equation is satisfied: 

/ W[''\x) = p[''\x,Y{x)) 



(12) 



Y{x) = V'{x)-W^i\x). 
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Notice that p[^\x, y) is a rational function of y, of degree N . If V' were a rational 
function of ; 
of degree n: 



function of x (finite n in Eqn 11 1, then p'f'\x, y) would be a rational function of x 



n N 

Pi'\^, y) = -9sTr^ + T.T. — ^ P^- 

" j=i i=i ^ " ' 

Thus the loop equation would be an algebraic equation, i.e. Y{x) would be an alge- 
braic function of x. Determining the rational function P^\x,y) (i.e. determining 
all the coefficients Pij) is possible but very tedious, and in fact, it is better to char- 
acterize an algebraic function Y{x) by its singularities and its periods. In general, 
one would find that w'^\x) = V'{x) — Y{x) has no singularity at the singularities 
of V , and the inverse function x{Y) has poles of residue gs at the eigenvalues of 
R. The genus of the algebraic function Y{x) and the periods § Ydx are related to 
the integration path C. 

In general, the relationship between the periods and C is quite complicated, 
but, for many applications to combinatorics, we are considering only formal matrix 
integrals, i.e formal perturbation with parameter i of a gaussian integral. In that 
case, the spectral curve Y{x) is always a genus curve (for the Hurwitz matrix 
integral, we prove it in Appendix|7|. This means that there exists a parametrization 
of Y{x) with a complex variable z: 

X = x{z) 
Y^y{z), 

where x and y are two analytical functions of z. The functions x{z) and y{z) are 
monovalued functions of z, but Y(x) — y{z) is multivalued, because there might 
exist several z such that x{z) — x. One of the determinations of Y{x) is called the 
physical sheet. 

The functions x and y are fully determined by their singularities. More precisely: 

• W^\x) = y'ix) — Y{x) is analytical in the physical sheet; it can have no 
singularity except at branchpoints. 

• From the definition of wf^ , we have that at x = oo in the physical sheet 

x\z) x\z) 

• As a consequence of Eqn |12[ x(^z) must have simple poles when z — > such 
that y(zi) = i.e. 

/ X as 

[z ~ z,} y'[Zi) 
In our case, for a potential of type [TT] this implies: 

^ x(z) = z + C-g.Eti 



Y{x) 



(13) 



^ {z-Zi)y'(zi 



^j=l (Z-Zj) x'(zj) 

where x{zj) — Xj, y{zi) — Ri. 

The above characterization of the spectral curve Y{x) is valid even if the potential 
V' is not rational, for instance if n — > oo. From now on, the function Y{x), or more 
precisely the pair of analytical functions 

S = z^ {x{z),y{z)) 
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is called the "spectral curve" of our matrix model. 



3.3. Topological recursion. We recall in this section the construction of the topo- 
logical recursion from For our purposes, we only deal with spectral curves of 
genus 0. Hence, we adapt the definitions of 0] to this case. 

A spectral curve 5 is a pair {x,y) of analytic functions on CPi. Let be the 
zeroes of dec, and assume they are simple. Then, locally at a^, y cx ^y{x — x{ai)). 
Let us denote z ^ z, the unique point corresponding to the other branch of the 
squareroot, such that x{z) = x{z). z is defined locally near the a^'s. 

A tower of fc-forms wjf\zi, . . . ,Zk) is constructed as follows. 

;(0) 



• Wl' = -ydx. 

• ^ is defined as the Bergmarj^ kernel: 

>Vf)(zi,Z2)=S(zi,Z2) = 



dzidz2 
(zi - Zlf 



• We define a recursion kernel 



K{z\z) 



{y{z) - y{z))dx{z) 



• For fc + 2g — 2 > 0, we define recursively the fc-forms VVjf\zi, . . . , Zk) by: 
Wi'\z,,z2,...,zk) = ^ Res if(zi,z)[w^^VH-^'^.^) 



K 



(14) + E <i+i(-'^)w.4/i(^'^\^) 

JCK, 0<h<g 

where ranges over (J, h) ^ (0, 0), (/, g). 

• = Tg are defined for 5 > 2 by 

^^-2^E.^-K(^w^) 

where $ is a primitive of ydx locally at the a^'s, i.e. — ydx. 

• The definition of J-i and J-q is more involved, and we refer the reader to 



(9) 



3.4. Main properties of the Wj 

• Symmetry Vfc,^, wl^\zi, . . . ,Zk) is symmetric in zi, . 



■ , Zk- 



• Invariance yk,g,2g + fc — 2 > 0, 
rational function of x. 



wlf'^ is unchanged if we add to y a 



Exchange invariance > 2 J-'g is unchanged if we exchange x and y. 



"^After Stefan Bergman (1895-1977), mathematician of Polish origin. 
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• Deformation If we perform an infinitesimal deformation (Sx, 5y) of the 
spectral curve, the W^.^^'s change. Let us introduce 

ri(z) = 5x{z)dy{z) — 5y{z)dx{z) . 

This form does not depend on the parametrization z, and can always be 
represented as /^/^^ B{z, z')A(z') for some path 7 and some meromorphic 
function A defined in its neighborhood (this data is called the dual of fi). 
Then 

(15) 5wi'\z^,...,zk)^ f A(z)>V^^i(z,zi,...,Zfc). 



• Limits wjf\s) is compatible with limits of curves. 

• Link to matrix models It was proved in [4J that, if our matrix model has 
a topological expansion property: 

i=l ^ ^' I c 9=0 

and 

00 

lnZ-^g2,-2^^^ 

3=0 

then, loop equations imply that: 

Wi'\zi, . . . , Zk)dx{z^) . . . dx{zk) = W['\zu. . . , zfc) 

where are computed using the spectral curve S = {x{z),y{z)) pre- 

sented in Section 3.2 (except for {k,g) = (1,0), (2,0), which receive simple 
additional contributions). Similarly, for g >2: 

F — T 

(Fq and Fi also receive simple corrections). 

Since the J-g(5)'s are invariant under transformations of S which leave |da; A dy\ 
unchanged, we call them symplectic invariants of S. 

In the following section we shall apply the general theory of to our Hurwitz 
matrix integral Eqn [8] 



4. Spectral curve of the Hurwitz matrix model 

In our case, the spectral curve of our matrix model must be determined order 
by order in powers of gs and t. Here, we only give the result. The proof is quite 
technical, and is deferred to Appendix [7] It relies on the fact that, after a suitable 
shift, to leading order in gg and t, we have a Gaussian matrix integral. In particular, 
this implies that the spectral curve has genus 0, i.e. it can be parametrized with 
a uniformization variable z € C. By composition with some homographic map, we 
can choose a parametrization for which x{z) ~ z when z 00 and x(0) — 0. 
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Since our potential V is of the form Eqn. [TT] with n oo, our spectral curve 
is of the form Eqn. |13[ Therefore we guess that the spectral curve must be of the 
form: 
(16) 

{x(z) = Z + Qs V . -I 7 ^ \ — 
^ ' ' ^1=1 {z—Zi)yi Ziyi 
viz) = -z + ln (z/t) +co + f-EZ^^f {Mz) - ^) 

where 



y{zi) ^ Ri = Invi 



and where 



is such that x{z) ' — has a finite limit when z — > 0. Cq is chosen such that 
V\x{z)) - y{z) - 0{l/z) at large 

1 ^ 1 

co = (A^-2)5.+.9.E— ■ 

The coefficient ci is 

_ 9s sr^ g, 
ci-y-2.^^/2M, 

(=1 

and one can check (this comes from the fact that the sum of all residues of V'{x) 
must vanish) that it is such that V'{x{z)) — y{z) ~ Nqs/z at large z i.e.: 

1 . x-^ 1 - Zi 



Each term Zi,yi,co,ci, is to be viewed as a power series in gs and in t. To the 
first few orders in g^ we have: 

N 



Co = gs 



ci = ^ + 0(5^ 
and notice that L{tvi) is a power series in t: 



i(t^.) = E 



m— 1 



The proof that 5(p; i) is the correct spectral curve for our problem is given in 
Appendix [7] It is obtained by computing the spectral curve order by order in the 
small gs and t expansion. 



The computation of wf,^^ 's in the topological recursion formula Eqn 



14 



involve 



taking residues at all zeroes of x'{z), i.e. involves symmetric rational functions of 



A MATRIX MODEL FOR SIMPLE HURWITZ NUMBERS, AND TOPOLOGICAL RECURSIONS 

the Zi's, and one can see that the coefficients of W^^'' are, order by order in gs 
and t, polynomials of pm — Qs X^i required for the computation of Hurwitz 

numbers indeed. 

4.1. Spectral curve at = 0. At = 0, the spectral curve reduces to: 

x{z) — z 

y{z) = —z + In (z/t) 

i.e. X = L{tey), where L is the Lambert function. Up to exchanging the roles of x 
and y, this is the Lambert curve 5Lambcrt appearing in Bouchard-Marino 



(17) 5o 



4.2. The symplectic invariants. So far, from the general theory of matrix models 
and from general properties of the topological recursion, we have proved that the 
generating function of simple Hurwitz numbers of genus g is the symplectic invariant 

(18) F,(p;.).X:*" E^ (2g-2+'; + %))! ^^-^-^^('^(P^^))' 



where 5(p; t) is the spectral curve of Eqn 16 (in fact we have proved it only for 
g > 2, and we consider that the cases g = and .g = 1 are easier). 

The computation of symplectic invariants J-g involves computing residues at 
the zeroes of x'{z), and there are 27V such zeroes, which makes the computation 
complicated. 

One may use the invariance properties of J-'g, under the exchange of x and y. 
Indeed, the zeroes of y'{z) are much simpler to compute, order by order in powers 
of (7s- At = 0, y'{z) has only one zero located at z = 1. 

Therefore, we introduce a new spectral curve satisfying 

Fgip;t)^Tg{S{p;t)), 
defined by exchanging x and y: 



S{p,gs;t) 



x{z) - -z + In (zjt) + Co+f-EZi /2i(^) 
y{z) = z + g, J2Zi i^=k)y^ + 5^ 



where now x{zi) — Invi and yi — x'{zi). 

4.3. The correlation forms. So far, we have introduced the Fg's as generating 
functions which encode the simple Hurwitz numbers Hg^^ by expansion on a proper 
basis of polynomials of the Pm's- Bouchard and Marifio [T define another generating 
function, namely the function of fc-variables: 



(2(7-2+H + fc)! '''''' 



where M^(x) = X^o-eSt Yii=i (un-normalized) symmetric monomials. 

It is easy to relate both. If we recall the combinatorial definition Eqn [18] 

we see that 
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We recall that Ri = Invi. 

The deformation property Eqn[T5]of symplectic invariants allows us to calculate 
their derivatives. When we perform an infinitesimal variation Vi Vi + 6vi, i.e. a 
variation Ri Ri + SRi on the spectral curve, we need to compute 

fli{z) = dx{z)dy{z) — Sy{z)dx{z). 

First notice that the form ydx is a meromorphic form (its singularities are poles) . 
It has simple poles of constant residues gs near z = Zj, i.e. locally near Zj we have 



y 



9s 



X — Rj' 



which implies that, locally near Zj, 

= ~9s ^ijSRj 



dx 



{x-R,Y 



0(1) 



Then, observe that the other poles of ydx are independent of Ri. For example, 
{V {x) — y)dx has a simple pole at oo, with residue Ng^ independent of Ri. There- 
fore, ri; has no residue, and thus no pole at oo. Similarly there is no pole at z = 0. 
The final result is that fi; is a meromorphic form with poles only at Zi : 

1 5Ridz 

— n,{z) = —. 

9s {z~-z,)^y. 



rii(z) can be written in term of the Bergman kernel B{z, z') 

1 



dz dz' 



as 



Q,i{z) ^ -gsSRi Res B{z,z') — 

z'^z, [z' - Zi) yi 

Then, the general theorems about the Tg's, tell us that 



5Tn = 



gJR, Res w['\z') , , ^ , 
z'^z, [z' - Zi) ; 



1 



gsSR, 



dx(zj) 



and more generally 



Swi'\z,,...,Zk)=gsSR, Res Wlf_>,{z,, . . . , z,, z') 

z'^Zi 

The result is that 

_ / Wi'\zi,...,zk) ' 



,{9) 



(z' - Zi)yi 



1 d^Tg 

g^dvi...dvk \dx{zi)---dx{zk)' 



Ri=x{zi) ,. . . ,Rk=x{zk) 



As a final step, we take the limit gs — > 0. The limit of the spectral curve iS(v; t) 
is simply 

x{z) = —z + In (z/t) 
y{z) = z, 

i.e. it is the Lambert spectral curve: y — L(te^). 



5, 



Lambert 
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In other words, we have proved that the function H'^^^vi, . . . , Vk) is the correla- 
tion form W^.^'' of the Lambert spectral curve iSLambort^ 







('^Lambert ) ; 


dx{zi) ■ ■ •dx(zfe) 



where x{zi) = Inw^ = Ri, i.e. Zi = L{tvi). 

This is precisely the Bouchard-Marino conjecture. 



5. Relationship with intersection numbers, Kontsevich integral and 

elsv formula 



Notice that the Lambert spectral curve 

^Lambert 



x{z) = —z + In (z/t) 
y{z) = z 



has only one branchpoint (solution of x'{z) = 0), given by z = 1. This is the reason 
why the Bouchard-Marino conjecture is so efficient to compute Hurwitz numbers. 

Since the topological recursion for computing the W^^^'s and ^g's, involves only 
the computation of residues at the branch point, we may perform a Taylor expansion 
near z — 1 : Let us define: 

y = i + C 

and 



-2(a; + 1 + lnt) 



We have, in the limit C 



■ 0: 



3 



4 



^ (-1)"'C^ 

m>2 



m 



and we invert that expansion y = 1 + ^ ^ 

y = 1 - 2^ 



£1 

36 



270 



e_ 

6 .6! 



which we write: 



m>l 



In other words, the wjf'^'s and Tg's of the Lambert curve 5Lambcrt, are the same 
;(<;) 



as the yVi s and JFg's of the following spectral curve Sk'- 



Sk 



x(0 = -l-\nt-\e 

2/(0-i-2e + E™>iW2e 



This spectral curve is exactly the Kontsevich spectral curve for times tmS (see 0). 
Here the t^'s satisfy the following recursion t2 — H — — \, and for m > 4: 



(19) 



m+l 



1=2 



We form the following series: 



y> (2to+ 1)! t2m+3 ^™ 

^ m! 2 ^ 

m— 1 
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and 



5(^) = _ln(l_/(;^))= 



m— 1 

we find to the first orders: 

z 8z^ 8z'^ 



Then, it was found in [S] that: 



do 



do-\ \-dk—3g—3-\-k j — 1 rniH \-mj—do,mi>0 

i=l i=l \i=l i=l I Mg.fc 

where M.g,k is the stable compact moduh space of Riemann surfaces of genus g with 
k marked points, and Kj is the j'^ Mumford's tautological class, and -ipi = ci{Ci) 
is the first Chern class of the cotangent bundle at the i*^ marked point. 

In other words, just by looking at the Lambert spectral curve, we see that there 
is a relationship between the generating function for Hurwitz numbers of genus g 
with a monodromy of length fc, and the generating function for intersection numbers 
of tautological classes on Mg^k- 

This type of relationship is completely natural and expected. The link coming 
from the ELSV formula [3J, relating Hurwitz numbers to Hodge integrals: 

^ (2.9-2 + %) + H)! '^Vr' [ AV(1) 

where Ag(i) = Y^l^^^t^ is the total Chern class of the Hodge bundle E over 
Mg,n- In fact, the change of variables taking the Hurwitz generating function to 
the generating function of all Hodge integrals with a single A factor is given by the 
Lambert curve itself [SIEIIII]. Furthermore, Mumford's formula fl\ ITS] 

00 / 21-2 \ 

ch(E) = 5 + E Kz-i + E (-1) V^^^'-^-M 

allows one to express Hodge integrals in terms of tp class intersections. The time 
parameters tm appearing in Eqn [19] are closely related to these topics. Mironov 
and Morozov |14j have previously considered similar constructions. 

An other natural question arising is the relationship between the Bouchard- 
Marino conjecture and the cut-and-join equation [8]. Since they are both recursive 
algorithms for computing Hurwitz numbers, it seems likely that they should be 
related. In fact, they are, appropriately interpreted, completely equivalent. This 
topic, as well as a more detailed discussion of the deformation of the Kontsevich 
model to the Hurwitz model, are deferred to future papers. 

6. Conclusion 

With the integral representation of V{N) characters, it is possible to express in 
general the partition function Z of Hurwitz numbers as a matrix model. In the 
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case of simple Hurwitz numbers, we obtain a 1-matrix model with external field, 
whose spectral curve is found by solving the master loop equation: 



5(p;t) 



2n 



xiz) = -z + In izjt) + CO + f - EZi ' 
Our main results are: 



f2niz) 



H^3\vi,...,Vk) 



dx(zi) • • • da;(z/j) 



(5, 



Lambert 



where 



x{zt) = Invi 
y{zi) = L{tVi) 



It provides an algorithm, namely the topological recursion of matrix models, 
to compute the Hg^^ by the residue formula of Paragraph 3.3 with only one 
branchpoint involved. As a matter of fact, this recursion relation between sim- 
ple Hurwitz numbers is understood to be equivalent to the Laplace transform of 
the cut-and-join equation with help of the elsv formula. Besides, Z for simple 
Hurwitz numbers is the time evolution of the kp r-function, a fact agreeing with 
its one-matrix model representation. We also see explicitly on the Lambert curve 
iSLambert = liiiigs^o '^(p; i) the relation between Z and the Kontsevich r-function. 

We hope that our matrix model- minded methods could help investigating double 
Hurwitz numbers (where Z is a Toda r-function) and further. 
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7. Appendix: proof of the spectral curve 

The proof works order by order in and t, and it relies on the fact that, to 
leading order, we have a Gaussian matrix integral. 

7.1. Shift of the matrix model. We start from: 

Nl A(v) 7h„(c) 

where the potential V{x) is: 

V{x) ^~+g,{N~^)x+ iHgs/t) + in)x~ g. In {T{-x/g.s)) + Ct, 
and Ct does not depend on x. Order by order in g^ we have the Stirling expansion: 

V'ix) = -x + ln (x/t) + - ^) + I - E igr- 

1=1 

We need to compute this matrix integral in the small g^ and t expansion (up to 
a constant factoring out of the integral) . 
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First let us perform a shift 

M = M-R 
where R = diag(^i, . . . , Rn) is such that: 

V'{Ri) = Ri. 

The equation V'{Ri) = Ri has several solutions, we choose the one which is a power 
series in Qs and t: 

Ri = L{tv,) - + + • • • = E^^^mK)- 

and we choose the determination of the Lambert function such that L{tVi) has a 
small t expansion L{tVi) = tvi + t^vf + • • • = Y^^^i m™~^{tvi)'^ /m\. 

Then we have: 



where we decompose V{M) 

V(M) = V2(M) + V>3(M) 
into a quadratic function of M 

V2(M) = ^Tr I -M^ + g^Y. IT^^^^^ 

y ~^ R - J3s R - jgs 

and V>3 contains all the higher degree terms (we will not need explicit expressions, 
however, the interested reader can derive them easily). 

Notice that V2{M) and V>3(M) have a small gs expansion (for instance approx- 
imate the sum in V2 by a Riemann integral). 

Order by order in we have: 



TV! A(v) 



If we rescale M = ^Jg^ A, we have 

/ d^e-^=(^) E yL(^^3(^/v^))-. 

Since V>3(A/y^) = 0{^/gl), we see that order by order in powers of gs, we may 
exchange the sum and integral. Therefore: 



- (-1)™ /• 

^0 "'^ -/w^CC) 



dMe-ii'"^^'^^ (V>3(M))'", 
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More generally, expectation values of polynomials Qp{M) are computed as formal 
power series, whose coefficients are polynomial moments of a gaussian integral: 



< Qp{M) > 



^m=0 m! JHn{C) 

In this form, we can use Wick's theorem. It shows that, if Qp{M) is any homoge- 
neous polynomial of total degree p in the entries of the matrix M , the expectation 
value (^Qp{M — R)^ is a power series in and t. It is expressed as a finite sum of 

connectecj^ fat-graphs. When we restrict the sum to the fatgraphs of genus g, we 
note it with a superscript We claim that, for our matrix model, this expectation 
value is 0(5?^^), and Oi^gl^-^). It can be seen as follows : fatgraphs contributing 
to the sum have one vertex of degree p, and v internal vertices of degree > 3 coming 
from (V>3)''. If we call / the number of faces, and e the number of edges in the 
fatgraph, its Euler characteristic is 2 — 2(7 = {v + \) — e + f . The power of gs 
coming from the gaussian integral is j > e — v (indeed, we have coming from the 

gaussian integral with e~'s^^^^^^\ and gj" accompanying (V>3)'', and in addition 
V2 and V>3 themselves have a gs expansion). So, we have: 

J >e-^; = 25-1 + / > 2.9-1. 

On the other hand, the number of half-edges is 2e = p + irii where rii is the 
number of internal vertices of degree i > 2, and we have v — n.i. So: 

EA . p p 

(2 - l)*^' +2-2' 

i 

In particular, let us show how to define the w[^'^''s, the topological expansion of 
the one-point correlation function. Consider : 

\ LV.T-R J x-Yl\ 
It is a double power series, whose coefficients are rational functions of x, and are 
polynomial gaussian expectation values of M. In its representation as a sum over 
fatgraphs, we collect those of genus g to define Tp^\x). Its coefficients are still 
rational functions of x. Since Tp^\x) = 0{g^^'^), we write 

Ti,'\^)= E 9lTl']{x), 
and thus we have, in the sense of formal power series of gg. 

00 2j 
p=0 J>2g-1 p=0 



We are in position to define W^^^ (x) as 



2j 



gi^-^wi^\x)^ V giy^T^hx). 



j=2g-l p=0 



■^Connected because the normalization factor Z ^ is included in the expectation value. 
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W^^"^ (x) is thus a formal power series in powers of gs and t, whose coefficients are 
rational functions of x. These definitions give a meaning to the equality between 
formal double power series: 

oo 

W^{x) = Y,gja-^wi'\x), 

9=0 

where Wi{x) is the resolvent : 

W^i(x)"/Tr ^ 



x-M 



In a similar manner, one can define (7^^"^+'^ wjf\xi, . . . , Xk) as the formal dou- 
ble power series computing the sum over (c for connected) fatgraphs of genus g 
arising in the correlation function: 

.,..)'™^-(nT,(-i^ ' 

By construction we have, in the sense of formal series: 

00 

W,{x,, ...,xk) =J29'''^^' Wi:'\x^, ...,xk). 

9=0 

To sum things up, the correlation functions wjf\xi, . . . , Xk) can be defined as 
formal power series in g^ and t, such that the coefficients are rational functions of 
the Xi's. It is defined by collecting the fatgraphs of genus g in the Wick theorem's 
expansion of gaussian integrals. 

The loop equations of gaussian matrix integrals are well known, and they imply 
that the wjf '' satisfy the topological recursion of [4]. 

For M^;^'''', or more precisely for Y{x) = V'{x) — Wi^\x), the loop equations read 
(see H): 

v'ix)-Y ^ uy'(^)-y'(^^) 1 '''' 



x-M y - R 
^ N 00 ^ ^ / / , X \ (0) 



-9^^' Y -K^^'^^^x~jgs Y~R, \\jgs-M 

A'^ 00 00 ^ ^ 

Y^ + -9^ E E E — ; T,,,, 

i=l j=Op=Q •'^ 



(20) 



- ( I I M ■ 



(0) 



d9s - R V jgs - R./ 

where as usual < . >^°-' means that we shift M = 11 + M, and keep only the genus 
zero fatgraphs in the gaussian expectation value. Notice that the sum over j is 
absolutely convergent. Moreover, to a given order in g^, the sum over p is finite. 

This equation is sufficient to determine Y{x) order by order in g^. To leading 
order we find Y{x) = —x + In {x/t) + 0{gs). To subleading orders, we recursively 
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have to determine a finite number of coefficients of the gs expansion of Tj,p;i . Those 
coefficients are completely determined by the condition that, to each order in g^, 
is a rational function of x with poles only a,t x — Ri, and in particular it 
must have no pole at a; = j^s, or at Y{x) — Ri. 



7.2. Asymptotic expansion of the spectral curve. To leading order in gs, we 
have V'{x) = —x + In (x/i), and: 

Y{x) = + In {x/t) + 0{gs) 

This leading order spectral curve has a rational uniformization: 

x{z) = z 

y{z) = -z + \n{z/t) 

Since the higher order gs corrections to V'{x) are all rational, the corrections to 
Pi^\x,y) are also rational functions of x and y. This implies that all corrections 
to Y(x) can be written with the uniformizing variable z, i.e. to all orders in gs the 
spectral curve is of genus 0. 

To all orders in gs, the equation Eqn [20] is algebraic of degree + 1 in the 
variable Y, this implies that the function x{z) is rational of degree N+ 1. It is easy 
to see that it has simple poles at Zi such that y{zi) = Ri, and one simple pole 
at oo. Up to a homographical change of variable z, we assume that x(0) ~ and 
^ z at large z, i.e.: 



1 1 



^ (a; - z.{)yi Ziyi 

Moreover, one sees directly from Eqn [20| that the residue of xAy is gs, i.e. yi = 
y' (zi). The function y(z) starts to leading order in gs as y{z) — — z+ln (z/i)+0(5s), 
and all the higher gs corrections are rational functions of z. Let x = oo and 
Xj = jgs, j <E N be the singularities of V{x). For each X^, let us choose Zj such 
that x(zj) — Xj (and we must choose the value of Zj which has a small t power 
series expansion). Since V'{x) has a simple pole of residue gs a± x — Xj, one sees 



from Eqn 20 that ydx has a simple pole of residue gs at z — Zj. Another way of 



saying this is to write: 



where Cq is a contour surrounding oo and all the Zj's. 



Order by order in gs, using the Stirling expansion of V'{x), this gives the spectral 
curve S{p; t) of Eqn 16 
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